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Abstract 

If neutrinos are to play a relevant cosmological role, they must be essentially 
degenerate. We study whether radiative corrections can or cannot be responsible 
for the small mass splittings, in agreement with all the available experimental 
data. We perform an exhaustive exploration of the bimaximal mixing scenario, 
finding that (i) the vacuum oscillations solution to the solar neutrino problem 
is always excluded; (ii) if the mass matrix is produced by a see-saw mechanism, 
there are large regions of the parameter space consistent with the large angle 
MSW solution, providing a natural origin for the Am^^^ <C Am^j^ hierarchy; 
(iii) the bimaximal structure becomes then stable under radiative corrections. 
We also provide analytical expressions for the mass splittings and mixing angles 
and present a particularly simple see-saw ansatz consistent with all observations. 
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1 Introduction 



Standard explanations of the observed atmospheric and solar neutrino anomalies IQ 
require neutrino oscillations between different species, which imply that neutrinos are 
massive, with mass-squared differences of at most 10~^ eV^. On the other hand, if 
neutrinos are to play an essential role in the large scale structure of the universe, the 
sum of their masses must be a few eV, and therefore they must be almost degenerate. 
This scenario has recently attracted much attention p[-||5|. In this paper we will analyze 
under which circumstances the "observed" mass differences arise naturally (or not), as 
a radiative effect, in agreement with all the available experimental data. 

Let us briefly review the current relevant experimental constraints on neutrino 
masses and mixing angles. Observations of atmospheric neutrinos require z/^ — v^- 
oscillations driven by a mass splitting and a mixing angle in the range 



On the other hand, there are three main explanations of the solar neutrino flux deficits, 
requiring oscillations of electron neutrinos into other species. The associated mass 
splittings and mixing angles depend on the type of solution: 
Small angle MSW (SAMSW) solution: 



5 X 10"^ eV^ < /\ml^ < 10"^ eV^ , 
sin^ 2e„t > 0.82 . 




3 X 10"*^ eV^ < Am2 



sol 



< 10"^ eV^ 



4 X 10"^ < sin^ 2e,oi < 1.3 x 10 



,-2 



(2) 



Large angle MSW (LAMSW) solution: 



10"^ eV^ < Am'l^i < 2 x 10^^ eV^ 
0.5 < sin^ 2esoi < 1. 



(3) 



Vacuum oscillations (VO) solution: 



5 X 10"^^ eV^ < Am^^j < 1.1 x 10"^° eV^ 
sm'^2esoi > 0.67. 



(4) 
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Let us remark the hierarchy of mass sphttings between the different species of neutri- 
nos, Aml^i ^ ^f^atJ which is apparent from eqs.(0, This hierarchy should be 
reproduced by any natural explanation of those splittings. 

As it has been shown in ref. the small angle MSW solution is unplausible in 
a scenario of nearly degenerate neutrinos, so we are left with the LAMSW and VO 
possibilities. 

An important point concerns the upper limit on sin^ 29soi in the LAMSW solution. 
According to ref. the absolute limit sin^ 26^30; = 1 is forbidden at 99.8%, but with no 
indication about the tolerable upper limit. In this sense a conservative upper bound 
sin^ 29soi < 0.99 can be adopted. However, recent combined analysis of data, including 
the day-night effect , indicate that even the sin^ 29soi = 1 possibility is allowed at 
99% for 2 X 10-5 eV^ < Am^^^ < 1.7 x 10"^ eV^ ahhough there may be problems 
to define the upper limit on sin^ 26*^0; in a precise sense ^ . For the moment we will 
not adopt an upper bound on sin^ 29 sol, but later on we will show the effect of such a 
bound on the results, which turns out to be dramatic. 

Finally, the non-observation of neutrinoless double /3-decay puts an upper bound 
on the ee element of the Majorana neutrino mass matrix, namely 

Mee <B = 0.2 eV. (5) 

Concerning the cosmological relevance of neutrinos, as mentioned above, we will as- 
sume Yl^rrij^. = 0(eV). In particular, we will take J^f^ui = 6 eV as a typical possibility, 
although as we will see, the results do not depend essentially on the particular value. 
It should be mentioned here that Tritium /3-decay experiments indicate m,^. < 2.5 eV 
for any mass eigenstate with a significant z/g component |[T0[| . 

Using standard notation, we define the effective mass term for the three light (left- 
handed) neutrinos in the flavour basis as 

C = -h^'^MyV + h.c. (6) 

where A^^ is the neutrino mass matrix. This is diagonalized in the usual way 

Mu = V*DV\ (7) 

where V is a. unitary 'CKM' matrix, relating the flavor eigenstates to the mass eigen- 



2 



states 



C2C3 C2S3 
-C1S3 - SiSaCae*'' C1C3 - SiSaSsC*'' 







S1C2 




C1C2 / 


^3/ 



where Sj and Cj denote sin 6i and cos di respectively. D may be written as 

/mie'^ \ 
D = \ m2e^<^' . (9) 




























It should be mentioned here that, for a given mass matrix, the 6i angles are not uniquely 
defined, unless one gives a criterion to order the mass eigenvectors z/j in eq.(^ (e.g. 
ml^ < < ml,J. Of course, the corresponding V matrices differ just in the ordering 
of the columns, and thus are physically equivalent. 
In this notation, constraint (|^) reads 

Mee = \m^, clcle'^ + m^, c^s^e^"^' + m^, si e^^^l < B . (10) 

As it has been put forward by Georgi and Glashow in ref. 0, a scenario of nearly 
degenerate neutrinos plausibly leads to ^^2 — 0. In that case, and for my = 2 eV, 
eq.(|TO|) yields sin^ 2^3 > 0.99, which, as discussed above, might be in conflict with the 



LAMSW solution to the solar neutrino anomaly. However, according to some flts, 62 
could be as large as 27^ or even larger 0. Therefore, although small values of 62 are 
clearly preferred, in some acceptable cases a non-negligible contribution of sin^ 62 in 
eq.(p!OD is enough to relax the above mentioned stringent bound on sin^2^3, and in 
fact we will see some examples of this feature later on. Our criterion throughout the 
paper will be to keep eq.(|^) [or the complete eq.([TO|)] as the neutrinoless double /3-decay 



constraint, without demanding any extra condition on 63, 62 [0. In any case, we will 
see that in all viable cases sin^ 2^3 > 0.99 and 62 is very small. 

Let us now discuss the strategy we have followed to analyze if the required mass 
splittings and mixing angles can or cannot arise in a natural way through radiative 
corrections. 

Following ref. (some of their arguments have been discussed above), the scenario 
of nearly degenerate neutrinos should be close to a bimaximal mixing, which constrains 
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the texture of the mass matrix 7Vl,y to be 



Mt 



( 
1 

1 



1 


1 \ 
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~2 
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~2 


2 y 



111 



where is a general mass scale. M.^ can be diagonalized by a \^ matrix 

-1 



/zi _ 

\/2 V2 
1 1 



1 



V 



121 



2 2 72/ 

leading to exactly degenerate neutrinos: D = rriy diag(— 1, 1, 1) and 62 = 0, sin^ 26^3 = 
sin^ 2^1 = 1. Let us remind that in this scenario only the LAMSW and VO solutions 
to the solar anomaly are acceptable. 

The nice aspect of A^;, suggests that it could be generated at some scale by inter- 
actions obeying appropriate continuous or discrete symmetries. This is an interesting 
issue |jl2|], which we will not address in this paper. On the other hand, in order to 
be realistic, the effective matrix Ai^, at low energy should be certainly close to Aib, 
but it must be slightly different in order to account for the mass splittings given in 
eqs. (|l|J^,|) . The main goal of this paper is to explore whether the appropriate split- 
tings (and mixing angles) can be generated or not through radiative corrections; more 
precisely, through the running of the renormalization group equations (RGEs) between 
the scale at which A^;^ is generated and low energy. The output of this analysis can be 
of three types: 

i) All the mass splittings and mixing angles obtained from the RG running are in 

agreement with all experimental limits and constraints. 

ii) Some (or all) mass splittings are much larger than the acceptable ranges. 

in) Some (or all) mass splittings are smaller than the acceptable ranges, and the rest 
is within. 

Case (i) is obviously fine. Case (ii) is disastrous. The only way-out would be 
an extremely artificial fine-tuning between the initial form of Ai^, and the effect of 
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the RG running. Hence we consider this possibihty unacceptable. Finally, case (in) 
is not fine, in the sense that the RGBs fail to explain the required modifications of 
Ail,. However, it leaves the door open to the possibility that other (not specified) 
effects could be responsible for them. In that case, the RGE would not spoil a fine 
pre-existing structure. Gonsequently, we consider this possibility as undecidable. We 
will see along the paper different scenarios corresponding to the three possibilities. 

Goncerning the mixing angles, it is worth stressing that, due to the two degenerate 
eigenvalues of A4b, Vb is not uniquely defined (H times any rotation on the plane of 
degenerate eigenvalues is equally fine). Hence, once the ambiguity is removed thanks 
to the small splittings coming from the RG running, the mixing angles may be very 
different from the desired ones. However, if those cases correspond to the previous (Hi) 
possibility, they are still of the "undecidable" type with respect to the mixing angles, 
since the modifications of Ai (of non-specified origin) needed to reproduce the correct 
mass splittings will also change dramatically the mixing angles. 

In section 2, we first examine the general case in which the neutrino masses arise 
from an effective operator, remnant from new physics entering at a scale A. In this 
framework, we take a low energy point of view, assuming a bimaximal-mixing mass 
structure at the scale A as an initial condition. In this section we do not consider 
possible perturbations of that initial condition coming from high-energy effects. We 
find this case to be of the undecidable type [possibility (Hi) above], except for the VO 
solution, which is excluded. 

In section 3 we go one step further and consider in detail a particularly well mo- 
tivated example of the previous case: the sec-saw scenario. We include here the high 
energy effects of the new degrees of freedom above the scale A (identified now with the 
mass of the right-handed neutrinos). We find regions of parameter space where the 
neutrino spectrum and mixing angles fall naturally in the pattern required to explain 
solar (LAMSW solution) and atmospheric neutrino anomalies, which we find remark- 
able. We complement the numerical results with compact analytical formulas which 
give a good description of them, and allow to understand the pattern of mass spht- 
tings and mixing angles obtained. We also present plausible textures for the neutrino 
Yukawa couplings leading to a good fit of the oscillation data. Finally we draw some 
conclusions. 
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2 M-jy as an effective operator 

In this section we consider the simplest possibihty that the Majorana mass matrix for 
the left-handed neutrinos, J^^, is generated at some high energy scale, A, by some 
unspecified mechanism (we allow A to vary from Mp to Mz)- Assuming that the only 
light fields below A are the Standard Model (SM) ones, the lowest dimension operator 
producing a mass term of this kind is uniquely given by [jl4| 

- -Kv^vHH + h.c. (13) 
4 

where k is a matricial coupling and H is the ordinary (neutral) Higgs. Obviously, 
= \k{H)'^. The effective couphng n runs with the scale below A, with a RGE 
given by IjlSl 



^ = [-?>gl + 2A + er/ + 2TrYtYe] - ^ [«:YtYe + {YIY^^k] , (14) 



where t = log/i, and g2, \,Yt,Ye are the SU{2) gauge coupling, the quartic Higgs 
coupling, the top Yukawa coupling and the matrix of Yukawa couplings for the charged 
leptons, respectively. Let us note that the RGE depends on A, and thus on the value 
of the Higgs mass, mu. We have taken mu = 150 GeV throughout the paper, but 
in any case the dependence is very small (it slightly affects the overall neutrino mass 
scale but not the relative splittings). 

In the scenario in which we are interested (almost degenerate neutrinos), the sim- 
plest assumption about the form of k is that the interactions responsible for it produce 
the bimaximal mixing texture of eq.(|ll]). Hence 

M,{K) = ]^K{K){Hf = M,. (15) 

Clearly, the last term of the RGE (|I^) will slightly perturb the initial form of Aiy, so 
we expect at low energy small mass splittings, and mixing angles different from the 
bimaximal case. 

In order to gain intuition on the final (low energy) form of M.ui and the correspond- 
ing mass splittings and mixing angles, it is convenient to neglect for a while all the 
charged lepton Yukawa couplings but Y^-. Then k maintains its form along the running. 
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Figure 1: Dependence of neutrino mass splittings at low energy {Amfj in eV ) with the cut-off scale 
A(GeV). 



except for the third row and column: 

/ 



MuifJ^o) oc 



1 



1 

72 



V2 
1 



-^(1 + ^) 



(16) 



where /iq is the low-energy scale, which can be identified with Mz and e has positive 
sign. Therefore, the mass eigenvalues are proportional to 1, —1 — e/2, 1 + 3e/2, and 
the corresponding splittings (adopting the convention m^^ < m^^ < ''^L) 



Ami2 



-Amo 



23 



13 



m,,e. 



(17) 



Clearly, these mass splittings are incompatible with the required hierarchy Am^^^ ^ 
Am^^, of eqs.(^ ^ We will discuss the size of these splittings shortly. Equation 
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([T6|) is also useful to get an approximate form of the V matrix responsible for the 
diagonalization oi Ai^ 

V ^ 
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!_ \ 


72 
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1" / 



which leads to mixing angles 
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sin2 2^i = — , sin2 2^. 
25 



K sm^2^.^^. 



, , (19) 

9' "25 ^ ' 

Clearly, these values are far away from the bimaximal mixing ones. In consequence, 
they are not acceptable. 

The previous failures of the scenario considered in this section in order to reproduce 
the mass splittings and the mixing angles indicate that we are in one of the two pos- 
sibilities (a) or (in) discussed in the Introduction. To go further, we need a numerical 
evaluation of the mass splittings, and thus of e. Solving the RGE ([T^) at lowest order, 
we simply find 



(20) 



So, from eq.(P!7D the mass splittings are typically of order 10~^ eV^. This is confirmed 
by the complete numerical evaluation of the RGE , which gives the mass splittings 
shown in Fig.l. 

The first conclusion is that for any A (even very close to yUo) the mass splittings 
are much larger than those required for the VO solution to the solar neutrino problem, 
Am^o; ~ 10^^° eV^. Therefore, the effect of the RGE for this scenario is disastrous in 
the sense discussed in the Introduction for the possibility (ii). In consequence the VO 
solution to the solar neutrino problem is excluded. 

For the LAMSW solution to the solar neutrino problem, things are a bit different. 
Clearly, the mass splittings obtained from the RGE analysis are within or below the 
required range, eq.(|^), for almost any value of the cut-off scale A. In addition, the 
mass splittings are clearly below the atmospheric range, eq.(|I|). Therefore, we are in 
the case (ii) discussed in the Introduction: The radiative corrections fail to provide an 
origin to the required mass splittings and mixing angles, but they will not destroy a 
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suitable initial modification (from unspecified origin) of Ai,y at the A scale. It is anyway 
remarkable that the radiative corrections place Am^ just in the right magnitude for 
the LAMSW solution. 

In the next section we will discuss a natural source for the initial modification 
of Ai,^, which leads naturally to completely satisfactory (atmospheric plus LAMSW) 
scenarios. 



3 Aijy from the see-saw mechanism 

A natural way of obtaining small neutrino masses is the so-called see-saw mechanism 



[jl6| in which the particle content of the Standard Model is enlarged by one additional 
neutrino field (not charged under the SM group) per generation, h'a,R (a = e, /i,r). 
The Lagrangian reads 

£ = -z/^m£,z/ + ^z/^Xz/J + h.c. (21) 

Here is a 3 x 3 Dirac mass matrix with magnitude determined by the electroweak 
breaking scale, m^) = Y^{H), where Y^, is the matrix of neutrino Yukawa couplings 
and {H) = 246/^2 GeV; and M is a 3 X 3 Major ana mass matrix which does not 
break the SM gauge symmetry. Thus, the overall scale of M., which we will denote by 
M, can be naturally many orders of magnitude higher than the electroweak scale. In 
that case, the low-energy effective theory, after integrating out the heavy i/^,/? fields 
[whose masses are 0(M)], is just the SM with left-handed neutrino masses given by 

= mlM-^mo = Y/M-^Y,{H)\ (22) 

suppressed with respect to the typical fermion masses by the inverse power of the large 
scale M. 

In this appealing framework, the degeneracy in neutrino masses can come about 
as a result of some symmetry (at some high energy, e.g. Mp) in the textures of A4 
and Yj, such that Aiu{Mp) = Aib- Starting from this very symmetric condition at Mp 
we make contact with the low energy neutrino mass matrix using the renormalization 
group. First we run Ai and Yj, from Mp down to M, where the right-handed neutrinos 
are decoupled. Below that scale and all the way down to low energy we run A^;^ as an 
effective operator, in the same way as in the previous section. We can also think of 
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the first stage of running between Mp and M as tlie fiigli energy effects providing an 
starting form of Ai^, for tlie second stage of running, i.e. the only one considered in 
the previous section (the role of A is played by M). This form, being slightly different 
from A^ft, may rescue the previous "undecidable" cases. 

In pursuing this idea, it is natural to assume [0 that the structure leading to 
the relation Ai,y{Mp) = Ai^ occurs either in the Dirac matrix rriD [while the Majorana 
matrix Ai is simply diagonal with equal eigenvalues (up to a sign)0] or in the Majorana 
matrix [with a simple diagonal Dirac matrix (again with eigenvalues equal up to a 
sign)]. We cannot expect a conspiracy between both matrices (which come from totally 
different physics) leading to A4iy{Mp) = M.^. 

3.1 Textures of M. and Yy leading to bimaximal mixing 

In the analysis of what structures for M. or 'Yy lead to M.y{Mp^ = Ai^ for simplicity 
we only consider real matrices (thus avoiding potential problems with CP violation). 

If we assume that is proportional to the identity and all the structure arises 
from the Majorana mass matrix, then we simply need to impose Ai{Mp) oc Aib (note 
that Ai'i^^ oc Aib )• The alternative case is less trivial but equally simple. If Ai is of 
the form0 




M = M 10, (23) 



then it is easy to see that the most general form of Yy satisfying Yy'^Ai ^Yy oc Ait is 

= YyBV,^ . (24) 
Here Yy is the overall magnitude of Yy and i? is a combination of two 'boosts' 



cosh a sinh a 
B= \ 1 

sinh a cosh a 



( cosh 6 sinh 6 
sinh 6 cosh 6 I , (25) 



V 1 

with two free parameters a, h. Actually, one could also take = YyRBVf^ , where R is 
a rotation in the (/i, r) plane, but such rotations can be absorbed into a change of the 

^Actually, the only condition we are imposing is that the eigenvalues are equal (up to a sign). A 
suitable transformation will diagonalize A4 to the form we assume. 

■^If A4 is taken proportional to the identity, the Yukawa matrix Yy must be chosen complex. Our 
choice of A4 avoids that complication and is physically equivalent. 
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{j^Rja basis, Rvr, with no modification in M. and thus give physically equivalent 

results. 

Let us notice that the former case, where all the structure is in the M. matrix, is 
equivalent to the latter case (all the structure in Y^^) if we set a = 6 = 0. To see this, 
note that a redefinition of the fields as vr — > Vi^vr, would leave M. oc diag(— 1, 1, 1) 
and Y,j = YyV^ . In consequence, it is enough to study the case where all the structure 
is in Yj,, given by eq.(|l). 



3.2 Running M.y from Mp to low energy 

From Mp to M the evolution of the relevant matrices is governed by the following 



renormalization group equations [ p!7[] : 
dY,, 1 _ 



dt 



IGtt 



2 ^ 



9 9 3 Q 
4^2 + l9i 



T I, 



-{y tY 



i e ^ e 



(26) 



C/Ye 

dt 



IGtt 



2 Ye 



+ '-^91 - t) Is + ^ (YjY^ 



i e ^ e 



where 



and 



T = Tr(3Y{jYu + SYJ^Yd + YIY, + Y+Ye), 



dM 



M{Y,Yjf + Y,YjM 



(27) 



(2^ 



(29) 



dt 167r2 

(not yet given in the literature). Here g2 and gi are the SU{2)l and U{1)y gauge 
coupling constants, and Yu,D,e are the Yukawa matrices for up quarks, down quarks 
and charged leptons. 

At M, vr decouple, and Ye must be diagonalized to redefine the flavour basis of 
leptons [note that the last term in ( pT]) produces non-diagonal contributions to Ye] 
affecting the form of the Yj, matrix. Then the effective mass matrix for the light 
neutrinos is My ~ Yj'M'^Yy{HY. 

From M to M^, the effective mass matrix A^,^ is run down in energy exactly as 
described in section 2. 

The renormalization group equations are integrated with the following boundary 
conditions: M. and Yy are chosen at Mp so as to satisfy 



My{Mp) = Mb, 



(30) 
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with the overall magnitude of fixed, for a given value of the Majorana mass M, 
by the requirement mi, ~ 0(eV). The boundary conditions for the other Yukawa 
couplings are also fixed at the low energy side to give the observed fermion masses. 
The free parameters are therefore M, a and h. 

3.3 Limits on the parameter space 

We discuss here the limits on the parameter space of our study, which is expanded by 
M, a and h. 



The parameters a and 6, which define the texture of the Y^, matrix through eqs. 



25), can be in principle any real numbers. However, it is apparent from eqs. (24, 25) 
that if a or 6 are large, the entries of Y,^ are extremely fine-tuned. Notice that the 
relative difference between cosh a and sinha factors is ~ 2e~^'"' (and the analogue for 
cosh h and sinh h). Therefore, the Yj, matrix is fine-tuned as ~ 

2g-2{|a| + |b|) pg^^^g 

m one. 

In particular, if \a\ or |6| are > 1.5 the matrix elements are fine-tuned at least in a 
10%. In consequence, we will demand 

|a|,|6|<1.5. (31) 

Let us remark that the previous limits are based on a criterion of naturality for the 
Yj, matrix. However, let us mention that if we relax these limits, the final results (to 
be presented in the next subsection) are basically unchanged, since the allowed regions 
for a, h remain in the "natural" region ( pi] ) in most of the cases. 

Concerning the remaining parameter, M, there is an upper bound on it coming 
from the fact that for large values of M, the neutrino Yukawa couplings, Yj,, develop 
Landau poles below Mp, spoiling the perturbativity of the theory |T8|. This occurs for 
M[Mp) ~ 4.3 X 10^^ GeV. Actually, there is an additional effect, namely the closing of 
the allowed Higgs window, which occurs approximately for the same value of M{Mp). 
Consequently, this sets the upper bound on M. It is interesting to note that this effect 
also restricts the values of a, b: for a given M, the larger a, b, the larger the entries of 
Yj,, and thus the lower the scale at which the Landau pole appears. In general, this 
restriction on a, b is less severe than eq.(pl|). 

Regarding the lower bound on M, there is no physical criterion for it, as the actual 
origin of the right-handed Majorana matrix is unknown. Since M can be written in 
terms of m,^ and Y^, (roughly speaking M ~ {H)"^ / niy) , we can adopt the sensible 
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Figure 2: Left plot: contours of ATOf2/eV^ in the (6, a) plane from loss than 10 ^ (black area), 
through 2 x 10~^ and 10~^ (hues) to more than 2 x 10~^ (grey). Right plot: same for Am^g/eV^, 
from 5 X 10"^ (black) to IQ-^ (grey). The Majorana mass is 8 x 10^ GeV. 

criterion that Y^, is at least as large as the smallest Yukawa coupling so far known, 
i.e. the electron one. This precisely corresponds to M ~ 100 GeV, below which is 
unplausible to descend. In consequence, our limits for M are 

lO^GeV <M < 4.3 X lO^^GeV. (32) 

On the other hand, as we will see in the next subsection, there are no physically viable 
scenarios unless 

M > lO^GeV, (33) 
which sets an operating lower limit on M. 

3.4 Results 

Figures 2 to 6 present our results for the mass splittings and mixing angles at low 
energy, after numerical integration of the RGBs from Mp to low energy as described in 
subsection 3.2 (we follow the convention ml^ < ml^ < ml^ in all the figures). We have 
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Figure 3: Left plot: Contours of sin^ 29^ in the (6, a) plane. The grey area marks the sin'^ 29^ > 0.64 
region. The line singled-out corresponds to sin^ 262 = 0.36. Right plot: Contours of sin^ 29i in the 
(&, a) plane. In the grey area sin^ 29i is smaller than 0.82, and the line corresponds to sin^ 29i = 0.9. 
The Majorana mass is 8 x 10^ GeV. 

chosen M = 8 x 10*^ GcV and m,y ~ 2 eV as a typical example; the dependence of the 
results with M and is discussed later on. 

Figure 2, left plot, shows contour hues of constant Am^2 (the squared mass dif- 
ference between the lightest neutrinos) in the plane {b, a) . The black (grey) region is 
excluded because there Ami2 < 10~^eV^ (Amfg > 2 x 10~^eV^), which is too small 
(large) to account for the oscillations of solar neutrinos (LAMSW solution). The white 
area is thus the allowed region. The hnes in it correspond to Amfa = 10~^ eV^ and 
2 X 10-5 eV^ 

Figure 2, right plot, gives contour lines of constant Am^g. The black (grey) region 
is excluded because there Am23 < 5 x 10^'' cV^ {Am^-!, > 10~^eV^), which is too small 
(large) to account for the oscillations of atmospheric neutrinos. Again, the white area 
is allowed. (The black area corresponds in fact to the "undccidablc" case discussed 
in the Introduction: it might be rescatable by unspecified extra effects.) We do not 
plot Amfg because it can always be inferred from Amig and Ami2- Moreover, in the 
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Figure 4: Left plot: Same as figure 3 for sin^ 263. The grey area corresponds to values above 0.99. 
The curve gives sin^ 263 — 0.95. Right plot: The grey area corresponds to cos 2^3 < 

interesting case, Amfg ^ ^^23 5 has Am^g ~ Am^g. 

The intersection of the white areas in both plots is non-zero and would give the 
allowed area concerning mass splittings. It is always the case that the area surrounding 
the origin is excluded. There, the mass differences are always of the same order, and 
follow the same pattern discussed in subsection 2 (Amgg = 2Am^2)- case we 

conclude that, away from the origin, there is a non negligible region of parameter space 
where it is natural to have Am^g ^ Amf2 and in accordance with the values required 
to explain the solar and atmospheric neutrino anomalies. In the following subsection 
we explain the origin of this hierarchy of mass differences. 

Next, we need to ensure that the mixing angles are the correct ones to give a 
good fit to atmospheric and solar neutrino data (as summarized by the ranges given 
in the Introduction). Figure 3, left plot, gives contours of constant sin^ 202 (one of the 
mixing angles relevant for atmospheric neutrino oscillations). The grey (white) area 
has sin^26'2 larger (smaller) than 0.64 and is disfavored (favored) by the data (SK + 
CHOOZ) at 99% CL. according to the most recent analysis (last paper of ref. 0]). 
The line singled-out corresponds to sin^26'2 = 0.36 (maximum allowed value at 90% 
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Figure 5: Region (two disconnected parts) in the {b, a) parameter space for M = 8 x 10^ GeV where 
all mass splittings and mixing angles satisfy experimental constraints. (See text for qualifications). 

C.L. according to the same reference). 

Figure 3, right plot, shows contours of constant sin^ 29i (the other mixing angle rel- 
evant for atmospheric neutrinos). The grey (white) area corresponds to sin^ 29i smaller 
(larger) than 0.82, and is thus disallowed (allowed). The additional line included has 
sin2 2^i = 0.9. 

Finally, figure 4, left plot, presents contours of constant sin^ 29^ which is relevant 
for oscillations of solar neutrinos. The grey (white) region has sin^ 29^ larger (smaller) 
than 0.99. If one is willing to interpret the existing data as impliying an upper bound 
of 0.99 on sin^ 26*3, then the grey region would be excluded. The plotted curve gives 
sin^ 2^^3 = 0.95 (sin^ 26'3 > 0.8 in all the region shown). Figure 4, right plot, shows 
the region of the parameter space accomplishing the resonance condition (cos 2^3 > 0), 
which is required for an efficient MSW solution of the solar anomaly (see however the 
first paper of ref. for caveats on this issue). 

The region of parameter space where all constraints on mixing angles and mass 
splittings are satisfied is given by the intersection of all white areas in figures 2, 3 and 
4 (right plot). If sin^2^3 < 0.99 is imposed, then that intersection region, including 
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Figure 6: 

Same as figure 5 for different values of tlie Majorana mass: Upper: 10^^ GeV; Lower left: 10^° GeV; 
Lower right: 10^ GeV. 
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now figure 4 (left plot), is empty and no allowed region remains. It should be noticed 
that this fact does not come from an incompatibility between the previous constraint 
and the sin^ 2^3 > 0.99 obtained from neutrinoless double /5-decay limits, eq. (0), in 
the 62 = approximation. If this were the case, it could be easily solved by decreasing 
the overall size of the neutrino masses, m^,, in eq . ([T0|) , and this is not the case. Indeed, 
eq.(|lOD is satisfied in nearly all the parameter space. Even where sin^ 29^ < 0.99, this is 
still true thanks to the contribution of 62. What actually forbids the whole parameter 
space if sin^ 26*3 < 0.99 is imposed, is the incompatibility between acceptable 9i, 62 and 
6'3 angles to fit simultaneously all the neutrino oscillation data, as can be seen from the 
figures. This fact remains when niy is decreased. In fact, the effect of decreasing nii, is 
essentially an amplification of the figures shown here, which comes from the fact that 
for a given Majorana mass, the neutrino Yukawa couplings become smaller (the effect 
is similar to decrease M, which is discussed below). 

If the sin^ 2^3 < 0.99 condition is relaxed (as discussed in the Introduction), then 
the allowed region is given by the two islands in figure 5, which is non- negligible. 

If we now vary M , this allowed region will move in parameter space as indicated 
in figure 6, where we show the allowed regions for a sequence of Majorana masses that 
range from 10^ to 10^^ GeV. As is apparent from the figure, lowering M has the effect 
of enlarging the allowed region which flies away from the origin, leaving at some point 
the region of naturalness for a and h. At M = 10^ GeV there is no allowed region inside 
the natural range for (a, h). Conversely, increasing M reduces the allowed region, which 
gets closer to the origin (at M = 10^^ GeV the allowed region becomes extinct). Let 
us remark again that in the allowed region one fits the observed atmospheric and solar 
anomalies, while part of the disallowed region corresponds in fact to the undecidable 
case (more precisely the region near the a = b = origin), in which some additional 
physics could be invocated to explain the same data. In contrast, the region away from 
the origin becomes excluded (the mass splittings are too large, even for atmospheric 
neutrinos). Let us also notice that if the sin^ 2^3 < 0.99 condition is imposed, the 
whole parameter space becomes disallowed for any value of M. 

We also find that, whenever there is a hierarchy in the mass splittings, the two 
lightest eigenvalues have opposite signs. This is just what is needed to have a cancel- 
lation occurring in the neutrinoless double /5-decay constraint ([T0|). This constraint is 
satisfied in almost the whole parameter space for any M. 
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This and other features of our results are discussed further in the next subsection. 



3.5 Analytical understanding of results 

It is simple and very illuminating to derive analytical approximations for the results pre- 



sented in the previous subsection. The renormalization group equations (p!4| , p6| , p7| , p9D 
we integrated numerically can also be integrated analytically in the approximation of 
constant right hand side. In this approximation (which works very well for our anal- 
ysis), the effective neutrino mass matrix at low-energy is simply Aif, plus some small 
perturbation. It is straightforward to obtain how the degenerate eigenvalues oi A4b 
get split by this perturbation. Neglecting the Y^,, Yukawa couphngs, we get the 
following analytical expressions: 



[-1 + {2clcl - l)e, - 2e 



1 + 3e^ - clcle^ ±\[^t + {clct - C2a)tu\ + 



1/2 



where = cosh a, S2a = sinh 2a, etc, and 

2 



1? 



1287r2 



1 ^^oi M^' 



34) 



(35) 



-^log^. 



(36) 



(The labelling of mass eigenvalues in eq.(^^ may not always correspond to the conven- 
tional order m^^ < m^^ < m1^.) It can be checked that, for a = 6 = 0, the Y,^ couplings 
play no role in the mass splittings. This is not surprising since, as was mentioned in 
subsection 3.1, this case is equivalent to having all the structure in A^, while Yj, is 
proportional to the identity. Then, it can be seen from the RGEs that all the non- 
universal modifications on A^j^ come from the Yg matrix, and has a form similar to 



the one found in section 2 [see eq.([l6|)]. So this scenario gives similar (not satisfactory) 
results to those found in that section. 

As ^ (which occurs as soon as Fj, > F,-, i.e for M > 10^ GeV), a further 
expansion in powers of er/e^ of the square root is possible in most of the parameter 
space (except where the coefficient of inside that square-root becomes very small). 
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The mass eigenvalues then read 



[-1 + {2clcl - l)e, - 2e.] , 



1 - {2clcl - l)e, + 2 




(37) 



■1/3 — 




These expressions show clearly the origin of the neutrino mass splittings. The splitting 
between the first two neutrinos is controlled by the small parameter e-r-, proportional 
to the squared Yukawa couplings of the charged leptons, and is insensitive to the 
parameter (proportional to the square of the larger Yukawa coupling Yy) which is 
responsible for the mass difference of the third neutrino. 

From the previous expressions we can extract the following remarkable conclusions. 
If the neutrino Yukawa couplings are sizeable (i.e. bigger than 1^), we will automati- 
cally obtain a hierarchy of mass splittings Am^g ^ Am23 ~ Am^g. This is exactly what 
is needed for a simultaneous solution of the atmospheric and solar neutrino anomalies, 
and thus represents a natural mechanism for the ^fn^oi ^ ^''^at hierarchy. Thus, if the 
mixing angle 62 is small, as it turns out to be in most of the parameter space, Am^2 is 
to be correctly identified with Am^^^ and Am^g with Am^^. The two mass eigenvalues 
which are more degenerate correspond to the lighter states, i.e. m^^ ~ m^^ < ""^^3- 
Moreover, m^-^ and 171,^^ have opposite signs in the diagonalized mass matrix, which 
is exactly what is needed to fulfill the neutrinoless double /3-decay condition, eq. ([TI1|) . 
All these nice features are illustrated by the explicit results presented in the previous 
subsection. 

Concerning the mixing angles, it is straightforward to check that, working in the 
Yi, > Yr approximation, the eigenstates of the perturbed matrix, corresponding to 
the previous eigenvalues, are 



where Vi are the eigenstates corresponding to the bimaximal mixing matrix, VJ, [see 
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and a, (3 are given by 
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(39) 



CaSfc + 0(e^/ei.), /3 = Sa + 0(e^/e^). 



(40) 



The Vl vectors define the new "CKM" matrix V' from which the mixing angles are 
extracted. Clearly, if just one of the two (a, h) parameters is vanishing, then V = Vb, i.e. 
exactly the bimaximal mixing case. Also, whenever Ca, Cb are sizeable (i.e. away from 
a = b = 0), \a\ ^ and thus we are close to the bimaximal case. Therefore, it is not 
surprising that in most of the parameter space shown in the previous section, this was 
in fact the case. This is remarkable, because it gives a natural origin for the bimaximal 
mixing, which was not guaranteed a priori due to the ambiguity in the diagonalization 
of the initial Ai^{Mp) = Aib niatrix, as was explained in the Introduction. 

On the other hand, the MSW condition cos 26*3 > (written using the conventional 
order m^^ < m^^) will be clearly satisfied as long as Vi corresponds to the lightest mass 
eigenvalue. In other words, this condition requires that the negative mass eigenvalue, 
see eqs.(|3^, [37| ), corresponds to the lightest neutrino. 



3.6 Examples of acceptable ansatze 

At a generic point in the allowed regions we have found, the form of the matrix Y^(Mp) 
would look rather ad-hoc: different elements in that matrix seem to conspire to give 
the correct neutrino mass texture. However, there are particular cases in which this 
matrix has a plausible structure. We give an example of such a texture for Yy{Mp) 
which for the case M ~ 8 x 10^ GeV studied in previous sections, would give mass 
splittings and mixing angles in agreement with observation: 
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(41) 



It corresponds to a = and b = sinh ^(3/4) ~ 0.69, value which falls in the allowed 
region plotted in figure 5. More precisely, the mass splittings are 

Am^2 - 2 X 10"^ eV^ Amj^ ~ 1 x 10~^eV^ Aml^ ~ 1 x 10"^ eV^ (42) 

and the mixing angles 

sin^ 202 = 0.04560, sin^ 20^ = 0.95405, sin^ 20^ = 0.99986. (43) 

Concerning the resonance condition for the MSW mechanism (see Fig. 4), this ansatz 
lies precisely at the border of the allowed area. 

Another examples of working ansatze can be obtained. For instance, the following 
ansatz (corresponding to a = cosh~"'^(-\/5/2) ^ 0.48, b = log(vT0/2) ^ 1.15 ) 

4 y/2 

1 V5 



Y,(Mp) = 



2v^ 2 2 
1 ^ v/5 



(44) 



V 4v^ 2 / 

works correctly for a wider range of right-handed Majorana masses (e.g. for M = 10^ 
GeV and M = 10^° GeV, as can be seen from Figure 6). 

It would be interesting to explore the possibility of finding a symmetry that could 
be responsible for the form of these ansatze and to analyze their implications for future 
long baseline experiments 

4 Conclusions 

We have performed an exhaustive study of the possibility that radiative corrections 
are responsible for the small mass splittings in the (cosmologically relevant) scenario 
of nearly degenerate neutrinos. To do that, we assume that the initial form of the 
neutrino mass matrix (generated at high energy by unspecified interactions) has the 
bimaximal mixing form, and run down to low energy. We then examine the form of 
the low-energy neutrino mass matrix, checking its consistency with all the available 
experimental data, including atmospheric and solar neutrino anomalies. 

We find cases where the radiative corrections produce mass splittings that are i) 
just fine ii) too large or in) too small. The vacuum oscillations solution to the solar 
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neutrino problem always falls in the situation {ii) , and it is therefore excluded. On the 
contrary, if the initial bimaximal mass matrix is produced by a see-saw mechanism (a 
possibility that we analyze in a detailed way), there are large regions of the parameter 
space consistent with the large angle MSW solution, providing a natural origin for the 
^"^loi ^ ^"^atm hierarchy. Concerning the mixing angles, they are remarkably stable 
and close to the bimaximal mixing form (something that is not guaranteed a priory, 
due to an ambiguity in the diagonalization of the initial matrix). 

We have explained analytically the origin of these remarkable features, giving ex- 
plicit expressions for the mass splittings and the mixing angles. In addition, we have 
presented particularly simple see-saw ansatze consistent with all the observations. 

Finally, we have noted that the scenario is very sensitive to a possible upper bound 
on sin^ 26'3 (the angle responsible for the solar neutrino oscillations). An upper bound 
such as sin^ 26*3 < 0.99 would disallow completely the scenario of nearly degenerate 
neutrinos due to the incompatibility between acceptable mixing angles to fit simulta- 
neously all the neutrino oscillation data. 

Addendum 

Shortly after completion of this work, a paper by J. Ellis and S. Lola on the same 
subject appeared |2^. In it, the scenario of our section 2 is also studied and similar 



(negative) conclusions reached. However, as we show in our section 3, positive results 
are obtained when the general see-saw scenario as the mechanism responsible for the 
effective neutrino mass matrix is studied. 

Also, the treatment by these authors of the constraints on mixing angles from 
neutrinoless double /3-decay and LAMSW fits to solar neutrino data is more pessimistic 
than ours. 
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